Capacitance of the Double Layer Formed at the Metal/Ionic-Conductor Interface: 

How Large Can It Be? 
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The capacitance of the double layer formed at a metal/ionic-conductor interface can be remarkably 
large, so that the apparent width of the double layer is as small as 0.3 A. Mean-field theories fail 
to explain such large capacitance. We propose an alternate theory of the ionic double layer which 
allows for the binding of discrete ions to their image charges in the metal. We show that at small 
voltages the capacitance of the double layer is limited only by the weak dipole-dipole repulsion 
between bound ions, and is therefore very large. At large voltages the depletion of bound ions from 
one of the capacitor electrodes triggers a collapse of the capacitance to the mean-field value. 



The rising demand for compact forms of energy storage 
with high power output has resulted in increased interest 
in electrochemical capacitors (ECs) [l|, 0]. An EC is a 
pair of metal electrodes separated by an ionic conductor, 
such as an aqueous solution of ions, an ionic liquid 3j, 
a super- ionic crystal [3], or an ion-conducting glass 5|. 
ECs with extremely high area per unit volume ( "superca- 
pacitors") are already used for a number of technologies. 
In this paper our focus is not on large surface area, but 
on the deeper physical question of a maximum possible 
capacitance per unit area. 

In a conventional double-plate capacitor, where metal 
electrodes of area S are separated by an insulator of 
width d and dielectric constant e, the capacitance C — 
eS/Aird (in Gaussian units). In an EC, the intervening 
medium is actually a conductor with finite conductivity 
cr, but with blocking of both ionic and electronic cur- 
rent at the electrode interface. The relation C — eS/And 
is therefore valid only at sufficiently high frequencies 
uj 3> 47r<7 /e, when the bulk of the ionic medium behaves 
as an insulator. We concern ourself with the opposite 
limit uj <C 47T(7/e, where polarization of the ionic medium 
eliminates electric field in the bulk and the capacitance 
of the EC is determined by the formation of thin electro- 
static double layers (EDLs) at both electrodes. 

How large can the capacitance be for these double lay- 
ers? The commonly-accepted expression for the max- 
imum possible capacitance of an EDL goes back to 
Helmholtz Q , who assumed that the charge of the metal 
surface is compensated by a layer of counterions with di- 
ameter a residing on the surface of metal. The resulting 
"Helmholtz capacitance" is given by 



C H = eS/2ira. 



(1) 



For a double-plate capacitor, where the EDLs formed at 
both electrodes can be thought of as two equal capaci- 
tances connected in series, the maximum capacitance is 
= eS/Ana. For a = 2A and e = 5, as we use below, 
C H /2S = 22 ^F/cm 2 . 

A recent experiment [5|, however, has reported much 
larger values of the EDL capacitance in phosophosili- 
cate glasses placed between platinum electrodes. Capac- 
itance per unit area as large as 400 /xF/cm was mea- 
sured, corresponding to an effective capacitor thickness 



d* = eS/AirC in the range 0.2 - 0.7 A, much smaller than 
any ion radius. The glass was held at a temperature of 
573 K, at which only the smallest ions, Na + with diam- 
eter a = 2A, are mobile. The dielectric constant of the 
glass e « 5. 

Current theories of EDL capacitors, based on the 
mean-field approach, fail to explain such large capaci- 
tance. The most widely-used theory is that of Gouy, 
Chapman, and Stern (CCS) 0, which extends the 
Helmholz capacitor concept to allow for the thermal mo- 
tion of counterions. In this approach, neutralizing ionic 
charge is imagined as a stack of thin uniform layers placed 
parallel to the charged electrode, with the charge density 
of each layer dictated by the Poisson-Boltzmann equa- 
tion. Such theories naturally lead only to a larger effec- 
tive capacitor thickness, and therefore a smaller capaci- 
tance than the Helmholtz value. 

In this paper we propose an alternate theory to ex- 
plain the large differential capacitance of the EDL. We 
abandon the mean-field approximation and deal instead 
with discrete ion charges, which interact strongly with 
the metal surface in a way that is not captured by the 
mean-field approximation. For simplicity, this paper will 
focus on the case where the conductivity a is due to a sin- 
gle species of ions with charge e. Our conclusions may be 
applied to capacitors made with ion-conducting glasses or 
superionic crystals, where only the smallest positive ionic 
species (such as Na + or Li + ) is mobile. 

An ion adjacent to a metal electrode produces elec- 
tronic polarization of the metal surface, and it experi- 
ences an attraction to the resulting image charge. For 
ions of small radius, the image attraction is significantly 
larger than the thermal energy fc^T, so that ions form 
stable, compact ion-image dipoles at the metal surface. 
Repulsion between adjacent dipoles results in the forma- 
tion of a strongly-correlated liquid of dipoles along the 
surface of both electrodes (see Fig. [TJ. When a voltage 
V is applied between the electrodes, the build-up of an 
electronic charge ±Q on the electrodes drives ions to de- 
tach from the positive electrode and to bind to the nega- 
tive one. Below, we demonstrate that the resulting EDL 
capacitance dQ/dV can be significantly larger than the 
Helmholtz value C#/2, since the dipole-dipole repulsion 



2 



that resists ion transfer is relatively small. 

Consequently, our answer to the question in the ti- 
tle is that the capacitance per unit area can be much 
larger than the Helmholtz value Ch/2S. In other words, 
the effective effective capacitor thickness d* can be much 
smaller than the ion radius. Below we derive an expres- 
sion for d* that is reasonably close to experimental values. 
Note that we do not assume any Faradaic effects, as are 
employed in theories of so-called "pseudo-capacitance" 
0, to explain large apparent capacitance values. 

Our theory also explains another peculiarity of the ex- 
periment 0, namely the sharp drop of differential ca- 
pacitance at a certain critical voltage. We show that, in- 
deed, capacitance should achieve a maximum at a partic- 
ular nonzero voltage before collapsing to a much smaller 
value. Contrary to the mean- field theory of Ref. @, this 
maximum is not driven by excluded volume effects among 
bound ions. Rather, the maximum occurs far below 
the complete filling of an ionic layer at either electrode 
@, when the voltage difference induces the positively- 
charged electrode to lose all of its bound ions. 

Here we consider the case of a parallel-plate capacitor, 
where the intervening medium is modeled as a fixed neg- 
ative background with charge density —eN, upon which 
resides a neutralizing concentration of mobile positive 
ions with charge +e and bulk density N. If we imag- 
ine the conducting ions to be hard spheres with diameter 
a, then an ion up against the metal surface experiences 
an attractive potential energy of approximately 



-e 2 /2ea. 



(2) 



For T = 573 K, e = 5, and a = 2 A, we get \u tm \/k B T w 
15, so that such ions are bound strongly to the surface. 
At a given voltage V, some area densities ni and ni of 
ions bind to the positive and negative plates, respectively. 
The following arguments generally assume that Na 3 < 1. 
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FIG. 1: A capacitor consisting of parallel metal plates (lightly 
shaded) bounding an ionic conductor with mobile positive 
ions. The neutral region of the ionic conductor is heavily 
shaded, while negative depletion regions (xi and x%) are left 
white. 

The attachment of positive ions to the metal plates is 
associated with the formation of regions with negative 



net charge q\ and (72 near the anode and cathode, re- 
spectively. Each of these charges exactly cancels the net 
charge of the adjoining plate and its bound ions, so that 
there is no electric field in the bulk of the ionic conductor. 
This implies q\ — —en\S—Q and 02 = —en^S+Q, where 
Q is the electronic charge that moves through the voltage 
source. In other words, the dipoles at each metal-glass in- 
terface are effectively embraced by a capacitor composed 
of the charge q and its positive image — q in the metal. 
As in every plane capacitor, the charges q% and 92 are 
uniformly distributed along the plane. 

Since positive ions gain a large energy \ui m \ by adsorb- 
ing to the metal surface, at equilibrium there must be a 
correspondingly large potential difference Vi m = \u im \/e 
between the metal surface and the bulk of the ionic con- 
ductor, so that the chemical potential of ions is uni- 
form. With such a large internal potential difference 
Vi m ^> ksT/e at each electrode, the negative regions 
are strongly depleted of ions even at very small applied 
voltage V. The charges q\ and 92 therefore constitute 
depletion layers of width x\ and X2 which form at the 
anode and cathode, respectively; here it is assumed that 
x\,X2 ^> a. These layers are devoid of mobile ions and 
have a charge density equal to that of the negative back- 
ground, so that qi = —eNxi and qi = —eNx%. Thus 



eNxi — erii — Q/S, 
eNx2 — en.2 = —Q/S. 



(3) 
(4) 



The electrostatic energy associated with the formation of 
the depletion layers can be estimated as 

U dep = [(«i + Q/Sef + (na - Q/Sef] . (5) 

In addition to Udep, there is a positive energy asso- 
ciated with repulsion between bound ions on the metal 
surface. When their density is low enough that na 2 <C I, 
ions repel each other by a dipole-dipole interaction: the 
potential created by a bound ion and its image charge 
repels an adjacent ion. In this limit the repulsive inter- 
action between two adjacent dipoles is approximately 



u d d{n) 



22 3/2 
e a n 1 



/2c. 



(6) 



This repulsion results in the formation of a strongly- 
correlated liquid of ions on the electrode surface, rem- 
iniscent of a two-dimensional Wigner crystal. The total 
dipole energy at a given plate is anSudd(n), where a is a 
numerical coefficient which depends on the structure of 
the lattice of dipole positions. For a triangular lattice, 
a sa 4.4. Further calculations will use this value. 

We now assemble a full description of the total energy 
U associated with the bound charge densities n\ and ri2- 
Taking as our reference the case where n\ = 112 = 0, 

U = S(ni+n2)u irn +Ude P +aS [n\Udd{ni) + n2Udd(n 2 )]-QV. 

(7) 

Here, — QV represents the work done by the voltage 
source. The equilibrium values of Q(V), ni(V), and 
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ri2(V) are those which minimize U. 
dU/dQ = gives 



The condition 



47T 



ma 



n 2 a 2 ) 



2Qa 2 
eS 



[ri2 — n\)a i 



1 in i 



Eq. ([5} suggests that at sufficiently small voltages, Q ~ 
eS(n 2 — ni)/2. Putting this result for Q back into Eq. 
we get a total energy U(ni,n 2 ). The conditions 
dU/dni = and 9n 2 = produce (via their addition 
and subtraction) the following relations: 



^3 (ma 2 + n 2 a 2 ) 2 + ^ [(ma 2 ) 3 / 2 + {n 2 a 2 ) 3 ' 2 



(n 2 a 2 ) 3 / 2 - (n ia 2 ) 3 / 2 = 



2V 

5aV in 



= 1, 
(9) 

(10) 



One immediate consequence of Eqs. © and (ITU)l is that 
at zero voltage, when the net charge Q of the capacitor is 
zero, there is still a finite concentration no of ions bound 
to each plate. If ions are relatively sparse in the bulk, so 
that Na 3 <C 1, then the second term of Eq. ((9]) can be 
neglected and 



n a 2 ~ ^jNa 3 /At: < 1. 



(11) 



Eq. (fTTjl implies that bound ions are sufficiently distant 
that a dipole-dipole interaction between them is justified. 
It also ensures that the depletion layer widths x±,x 2 a, 
as assumed earlier. 

As the voltage is increased from zero, ions are driven 
away from the anode and attracted to the cathode, so 
that ni decreases and n 2 increases. By comparing Eqs. 
(O and ([6]), we see that the condition Na 3 <C 1 im- 
plies that the energy cost associated with increasing the 
total number of bound ions, and thereby causing the de- 
pletion layers to swell, is much larger than the dipole- 
dipole interaction energy. As a consequence, the total 
number of bound ions n\ + n 2 ~ 2no remains almost 
constant with increasing voltage. The electronic charge 
Q ~ eS(n 2 — ni)/2 that passes through the voltage source 
can therefore be thought of as the corresponding move- 
ment of image charges from one plate to another. At a 
particular voltage V c , however, the anode becomes com- 
pletely depleted of bound ions. This occurs when n\ = 
and n 2 ~ 2uq, so that Eq. (|10[) yields 



Vc ~ 2 \ it ) Vm - 



(12) 



Substitution of n\ + n 2 = 2ng into Eq. © verifies that 
Q ~ eS(n 2 - m)/2 for all \V\ < V c . At \V\ > V c , the 
number of bound ions on the non-depleted electrode may 
still increase, but only through the costly widening of the 
depletion layer. 

Knowing ri\ and n 2 as a function of voltage allows us 
to calculate the differential capacitance 



r(v x _ dQ ^ eS ( dn 2 dn\ 
1 ' ~ dV ~ ~2 \~dV ~ ~d~V 



(13) 



Taking the derivative of Eq. (fT0|) with respect to V, and 
using n% + n 2 — 2n , gives 

^ 4 / 4tt \ lfi eS Svr / 4vr \ 1/4 . A . 
C{0 ^T^{Na- 3 ) - = T^{Na^) Cff(14) 

for the capacitance at zero voltage. For A^a 3 « 1 we 
arrive at C(0) ^> Ch/2. Such large capacitance arises 
because at V < V c charging of the capacitor is limited 
only by the dipole-dipole repulsion energy, which is small 
since no a 2 <C 1. 

Note that at V — the capacitor is composed of two 
identical EDLs in series, so that the total capacitance 
C(0) is equal to half the capacitance of each. At higher 
voltages, the capacitance of the EDL near the positive 
plate increases strongly as this plate becomes depleted 
of ions and the corresponding dipole-dipole energy goes 
to zero. Thus the contribution of the positive plate to 
the total capacitance vanishes at V = V c . As a result, 
immediately prior to V = V c the capacitance achieves its 
maximum value 

c --iL(a^) 1/4 T-^- ^ 

At V > V c the capacitance collapses to a much smaller 
value. This "depleted capacitance" can be found through 
optimization of the total energy U under the condition 
ri i = 0, which yields 



Cdep{V) — \ 



a V 47r(V/y im + l) 



(16) 



This expression neglects the weak dipole-dipole interac- 
tion at the non-depleted negative plate. At V/Vi m ^> 1, 
the capacitance is dominated by the depletion layer next 
to the positive electrode, and therefore it approaches the 
standard value for capacitance of a depletion layer. 

Fig. [5] shows the capacitance and the density of bound 
ions as a function of voltage, as calculated by a numerical 
minimization of the total energy in Eq. (JT)). We have used 
Na 3 =0.1, following the estimate of Ref. 

The effective thickness d* min corresponding to C max is 



15a { Na 



32tt 



a w 0.49 (Na 3 ) 1/4 a. (17) 



Thus we arrive at a remarkable prediction: the effective 
capacitor thickness can be much smaller than the ion di- 
ameter a. As an example, an ion-conducting glass with 
mobile sodium atoms of diameter a — 2 A and density 
Na 3 = 0.1 can be used to make a capacitor with capaci- 
tance nearly five times larger than Ch/2. As mentioned 
before, this surprisingly high capacitance is a result of 
the weak dipole-dipole interaction between bound ions 
that comprise the double layer. Indeed, near the capac- 
itance maximum the filling factor on the negative plate 
n 2 a 2 ~ 2n a 2 < 1, so that it is incorrect to think of the 
EDL as a series of uniformly charged layers. This differ- 
ence represents an important change of paradigm, from a 
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FIG. 2: The area densities of bound ions (dashed lines, left 
vertical axis) and the capacitance (solid line, right vertical 
axis) as a function of applied voltage for Na 3 = 0.1. 



mean-field capacitor to a capacitor composed of discrete, 
correlated dipoles. 

We can compare our theory to the experiments of Ref. 
Q. Capacitance- voltage characteristics for three differ- 
ent phosphosilicate glasses are shown in Fig. [3] together 
our theoretical prediction (dashed line), using e = 5 and 
N = 10 22 cm~ 3 as estimated by the authors of Ref. 
The agreement with theory can be improved (solid line) 
if we assume that the image attraction of Eq. ^ uses a 
smaller dielectric constant e m = 2 [? ], while the macro- 
scopic dielectric constant e = 10. 

While our theory seems to provide a qualitative ex- 
planation of the data, the observed collapse is not as 
sharp as we predict. This may be due in part to the 
fact that our theory is valid only in the limit Na 3 <C 1, 
and the use of a marginally small Na 3 = 0.1 does not 
allow for a large separation between the length scales a, 

— 1/2 

n Q , and X\, x%. We have also ignored temperature 
and disorder effects. The importance of thermal effects 
can be estimated by calculating the mean-square ther- 
mal displacement (r 2 ) 1 / 2 = 6 of a given dipole from its 
potential energy minimum. For e = 5 and T = 573 K, 
we get relatively small displacements 5 ~ 0.3n -1 / 2 at 
n = no and 5 ~ 0.2n -1 / 2 at n = 2tiq. Thus it seems un- 
likely that thermal motion alone can be responsible for 



the substantial smearing of the capacitance collapse. 

One can imagine that at much higher temperatures, 
such that fcflT 3> Udd(no) but ksT <C Ui m , the thermal 
energy of bound ions becomes much larger than their 
mutual repulsion, and therefore dipoles on a given metal 
surface are better described as a two-dimensional ideal 
gas than as a Wigner crystal. In this limit the capaci- 
tance is determined by the change in the ideal gas entropy 
associated with transfer of ions from one electrode to the 
other. This description gives a zero-voltage capacitance 
C(0) = n a 2 (uim/kBT)(eS/a). For the crossover tem- 
perature T c to Eq. (fT4"]) we get T c = 15audd(no) I '4fce ~ 
3700 K. Therefore, at T = 573 K, deviations from our 
low-temperature theory are unlikely to be large. The 
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FIG. 3: The capacitance of an ion-conducting glass between 
metal plates. The dashed and solid lines are theoretical pre- 
dictions, calculated by a numerical minimization of the energy 
in Eq. (J7J). The dashed line uses e — 5 while the solid line as- 
sumes e = 10, e m = 2. Thin lines with dots show data for 
three different glasses from Ref. §|. 



disorder potential acting on mobile ions in the glass may 
be more important, and we leave this for the subject of 
a later publication. 
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